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ABSTRACT 

A computer  program  which  calculates  the  inviscid 
approximation  of  propeller/stem  boundary- layer  interaction 
is  given.  The  description  includes  a summary  of  the  inviscid 
theory  and  the  assumptions  used,  an  outline  of  the  numerical 
procedure  and  computation  method,  a numerical  example,  FOR- 
TRAN listing,  and  input  and  output  instructions.  This  pro- 
gram is  designed  to  compute  the  effective  velocity  profile 
from  the  measured  nominal  velocity  profile  for  propeller 
design. 


ADMINISTRATIVE  INFORMATION 

This  work  was  sponsored  by  the  Naval  Material  Command/Naval  Sea 
Systems  Command  under  the  Direct  Laboratory  Funding  Program  for  Hydro- 
dynamics of  Very  High-Speed  Submarines,  Element  Number  62543N.  The  work 
was  performed  under  internal  Work  Units  1S20-004  and  1500-200. 


INTRODUCTION 

A propeller  axisymmetric- stern  boundary- layer  computer  program  is 
documented  on  the  basis  of  the  theory  derived  in  Reference  1.  The  para- 
meters relevant  to  this  computer  program  are  briefly  discussed  below. 

Many  ship  propellers  have  a radius  which  is  about  the  same  order  of 
magnitude  as  the  boundary- layer  thickness  at  the  propeller  location. 

The  nominal  wake  distribution  is  defined  as  the  velocity  distribution  at 
the  propeller  plane  in  the  absence  of  a propeller.  For  the  model,  the 


1 Huang,  T.T. , H.T.  Wang,  N.  Santelli,  and  N.C.  Groves,  "Propeller/ 
Stem/Boundary-Layer  Interaction  on  Axisymmetric  Bodies:  Theory 

and  Experiment,"  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  Report  76-0113  (Dec  1976) . 
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nominal  velocity  distribution  is  usually  measured  by  a standard  wake  rake 
at  the  propeller  disc  in  the  absence  of  a propeller.  The  nominal  wake 
distribution  for  the  full-scale  ship  is  usually  computed  by  adding  a 
scale-effect  correction  to  the  measured  model's  nominal  wake  distribution. 

An  operating  propeller  produces  an  upstream  suction  which  results  in 
an  increase  of  velocity  in  the  stern  boundary- layer  region,  as  shown  in 
Figure  1.  A typical  stream  surface  moves  inward  from  radius  r to  r^ 
while  the  velocity  is  increased  from  the  nominal  axial  velocity  u^  to 
Up.  The  resultant  axial  velocity  u^  in  front  of  an  operating  propeller 
is  called  the  apparent  velocity  distribution.  The  velocity  distribution 
actually  experienced  by  the  propeller  blade  in  developing  thrust  and  torque 
is  called  the  effective  velocity  distribution.  The  effective  velocity, 
which  is  to  be  used  for  the  design  of  a wake-adapted  propeller,  is  the 
apparent  velocity  minus  the  circumferential-mean  propeller-induced  axial 

velocity,  u . The  propeller- induced  velocity  can  be  computed  from  the 

^ 2 
propeller  field-point  velocity  program  developed  by  Kerwin  and  Leopold, 

if  the  loading  and  geometry  of  the  propeller  are  known.  However,  some- 
times only  the  propeller  geometry  is  known;  then  an  interactive  procedure 
to  compute  the  loading  (using  a propeller  inverse  computer  program  such 
as  the  one  developed  by  Cummings^),  propeller- induced  velocity  (refer- 
ence 2),  and  the  present  program  are  required.  A detailed  description 
of  the  theory,  computation  method,  and  iteration  procedure  is  given  in 
Reference  1. 


2 Kerwin,  J.E.  and  R.  Leopold,  "A  Design  Theory  for  Subcavitating 
Propellers,"  Transactions  of  the  Society  of  Naval  Architects  and 
Marine  Engineers,  Vol.  72  (1964),  pp.  294-335. 

3 Cummings,  D.E.,  "Numerical  Prediction  of  Propeller  Characteristics,"' 
Journal  of  Ship  Research,  Vol.  17,  Part  3 (1973),  pp.  12-18. 
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The  input  of  the  present  computer  program  is  the  nominal  axial 
velocity  profiles  at  the  plane  of  propeller.  The  present  program 

then  computes  the  apparent  (u  /u  ) and  the  effective  (u  /V  ■ (u_-u  )/V  ) 

P X 6 S p 3 S 

axial  velocity  profiles  based  on  the  following  assumptions: 

(a)  The  flow  is  axisymmetric,  the  fluid  is  incompressible,  and 
upstream  of  the  propeller,  the  mean  circumferential  velocity  is  identi- 
cally equal  to  zero  with  and  without  the  propeller  in  operation. 

(b)  The  interaction  of  propeller  and  stem  boundary  layer  is  con- 
sidered to  be  inviscid  in  nature.  Thus,  propeller-induced  viscous  losses 
and  turbulent  Reynolds  stresses  are  neglected. 

(c)  The  conventional  boundary- layer  assumption  of 

is  assumed  to  be  valid  for  the  boundary  layer  in  the  absence  of  a pro- 
peller; and 

(d)  Upstream  of  the  propeller,  no  energy  is  added  to  the  fluid  by 
the  propeller,  and  the  propeller- induced  velocity  field  upstream  of  the 
propeller  is  assumed  to  be  irrotational . 

Based  on  the  above  assumptions,  two  equations  governing  the  pro- 
peller/stern  boundary- layer  interaction  have  been  derived  in  reference  1, 


i.e. , 


r«S.<lr=  , 


M“;(= 


Since  u^  and  Up  can  be  approximated  locally  by  linear  functions  of  r 
and  r , the  mass  flux  within  the  stream  annuli  given  by  4r  * 


can  be  integrated  from  Equation  (1) . Integrat- 


ing Equations  (1)  and  C2) , one  obtains  the  finite-difference  forms  of 
the  governing  equations. 
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This  program  computes  the  values  of  u at  various  values  of  r from  the 

P P 

input  values  of  u (r)  and  u (r)  by  a simple  iteration  procedure.  The 
numerical  scheme  developed  for  solving  Equations  (3)  and  (4)  was  found 
to  converge  rapidly. 


NUMERICAL  PROCEDURE 

If  the  loading  and  geometry  of  the  propeller  are  known,  the  lifting-  f 

_ 9 

surface  computer  program  developed  by  Kerwin  and  Leopold"  can  be  used  to 
calculate  the  values  of  u at  arbitrary  points  in  the  boundary  layer. 

The  following  procedure  is  used  to  compute  the  values  of  u^  from  Equa- 
tions (3)  and  (4)  at  various  radial  locations  r^,  given  calculated 
values  u (r  ) . 

a P 

(a)  First,  assume  that  u^,  (rpi)  = uxi(rp.)  = u^.Crpp  for  i=l,N 
(i  = 1 on  body  surface,  i = N outside  the  boundary  layer). 

(b)  Use  a least-squares  fitting  technique  to  obtain  a polynomial 

» 

representation  of  the  measured  values  of  u (r.)  for  i=l,N. 

Xi  1 

4 * 


(c)  Use  Equation  (3)  to  determine  values  of  u (r  ) , under  the 

Pi  Pi 

assumptions  that  the  mass  flux  in  the  presence  of  the  pro- 
peller, from  r to  r , is  equal  to  the  flux  in  the  absence 
Pi  Pi+1 

of  the  propeller,  from  r^^  to 

fd)  Improve  the  values  of  u fr  ) by  u (r.)  + u (r  ) + u_  (r  ) 

p.  p.  ^ X.  a.  p.  C.  p. 

and  use  Equation  (4)  to  obtain  u^  (r^  ) . Since  outside  the 

boundary  layer  u (r^  .)  = u (r^)  = V , Equation  (4)  yields 

u (r  ) = u (r  ) ♦ V a u (r  ) + u (r„) . The  correction 

Pn  Pn  Pn  ^ Pn  ’'n 

velocity  u.  (r  ) for  the  second  approximation  is  zero  outside 
^i  Pi 

the  boundary  layer,  u (r  >6)  = 0.  For  r <6,  Equation  (4) 

^i  Pi  Pi 

can  be  used  to  solve  for  u_  (r  > 6)  step  by  step  from  the  edge 

^i  Pi 

of  the  boundary  layer  where  u (r  >^)  = 0 inwards  down  to 


i Pi 


the  radius  of  the  body  surface. 


(e)  Kith  the  new  set  of  r (r  ) repeat  steps  (c)  and  (d)  until  the 

Pi  Pi 

distribution  of  u versus  r converges. 


EXAMPLE 

Figure  2 shows  (a)  the  measured  nominal  axial  velocity  profile, 

u^/V^,  in  the  absence  of  a propeller  (Afterbody  1 of  Reference  1), 

(b)  the  final  propeller  nondimensional  circulation  distribution  G(r/R  ) 

3 P 

calculated  by  the  propeller  inverse  program,  (c)  the  final  effective 

velocity  (Ug/V^)  and  apparent  (Up/V^)  axial  velocity  distributions  com- 
puted by  the  present  propeller/ stem  boundary- layer  interaction  program. 
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and  (d)  the  propeller- induced  circumferential-mean  axial  velocity,  u /V 
computed  by  the  propeller  field-point  program.^  The  final  propeller^ 
thrust  coefficient.  » T / (i/»  is  0.371,  and  the  advance 

ratio,  J=Vj/Hl^,  is  1.25  for  this  example.  As  discussed  in  Ref- 
erence 1,  the  computed  values  of  and  G have  essentially  converged 

to  their  final  values  after  three  iterations. 
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INPUT  INSTRUCTION 

INPUT  STATEMENTS 

The  input  statement  by  means  of  which  data  are  entered  into  the 
program  are  as  follows: 

READ  (5.1)  NI 

READ  (5.7)  (TITL(I),  I * 1.8) 

READ  (5.1)  N,  NPRL 

READ  (5.5)  (UX(D)  R(I),  UAP(I),  I = 1,N) 

The  corresponding  format  statements  are  as  follows: 

1 FORMAT  (2110) 

7 FORMAT  (8A10) 

5 FORMAT  (3F10.5) 


DEFINITION  OF  INPUT  VARIABLES 

; MI  Number  of  input  sets  or  cases  to  be  computed 

(: 

n 

i . TITL  Title 


^ N Total  number  of  points  from  hub  out  to  the  edge  of  the  boundary 

i’  layer  where  nominal  velocities  are  given 

s 

r 

NPRL  Total  number  of  points  from  hub  out  to  the  tip  of  propeller 

where  nominal  velocities  are  given;  NPRL  may  be  smaller  than  N 


UX(I) 


u^/V^,  nominal  axial  velocity/ship  speed 


Radius  (use  dimensionless  radius  r/R^,  R^  being  the 
propeller  radius 


UAP(l)  u /V  , propeller- induced  axial  velocity/ship  speed 

£L  S 


DEFINITION  OF  OUTPUT  VARIABLES 

The  output  variables  are  defined  in  the  order  that  they  appear 
in  output  listing. 


UX(I) 


UPD(I) 


Nominal  axial  velocity  ratio,  (input  variable) 


Radius  (input  variable) 


Effective  axial  velocity  ratio,  Ug  /Vi 


UP (I)  Apparent  axial  velocity  ratio,  u /V  , the  resulting  velocity 

P ^ 

with  propeller  in  operation 


UDF(I)  Velocity  ratio  different  with  and  without  propeller  in 

operation  (u^/V^)  - (u^^/V^) 

UAP(I)  variable) 

WR(I)*  Volume-mean  nominal  axial  velocity  ratio,  /y, 

WPR(I)  Volume-mean  effective  axial  velocity  ratio  ^ 

^ S 


WR(I)/WPR(I) 
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iM*  continue 

STOP 
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Figure  2 - Computed  Nondirr.ensional  Circulation  and  Effective  and  Apparent 
Axial  Velocity  Profiles  after  Four  Iterations  for  a Typical  Propeller 
Operated  in  a Measured  Nominal  A<ia(  Velocity  Profile  of  a Typical 
Axisymmetrical  Body 
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(2)  DEPARTMENTAL  REPORTS.  A SEMIFORMAL  SERIES.  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE,  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION. 

(3)  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
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